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Abstract 

The quark-meson coupling (QMC) model, which has been successfully used to describe 
the properties of both finite nuclei and infinite nuclear matter, is applied to a study of A 
hypernuclei. With the assumption that the (self-consistent) exchanged scalar, and vector, 
mesons couple only to the u and d quarks, a very weak spin-orbit force in the A-nucleus 
interaction is achieved automatically. This is a direct consequence of the SU(6) quark 
model wave function of the A used in the QMC model. Possible implications and exten¬ 
sions of the present investigation are also discussed. 
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In earlier work we addressed the question of whether quarks play an important role in 
finite nuclei a q This involved quantitative investigations of the properties of closed-shell 
nuclei from 16 0 to 208 Pb. These calculations were performed within the quark-meson coupling 
(QMC) model, first suggested by Guichon ||, where the interaction between nucleons involved 
the exchange of scalar (er) and vector (a; and p) mesons self-consistently coupled to the quarks 
within those nucleons. Within the model it has proven possible to successfully describe not only 
the properties of infinite nuclear matter Q, but also the properties of finite nuclei |l], Q. Blunden 
and Miller || , and Jin and Jennings |J have made similar studies of the QMC model and some 
phenomenological extensions. One of the most attractive features of the QMC model is that 
it does not involve much in the way of additional complications to Quantum Hadrodynamics 
(QHD) [jj. Furthermore, it produces a reasonable value for the nuclear incompressibility [4|]. 
Recently, with an extended version of the QMC model which took account of the quark structure 
of the to and p mesons as well, we studied the properties of finite nuclei in a unified manner, as 
well as the hadron mass changes in finite nuclei 0 . 

Here we develop the model further to study the properties of A hypernuclei. One of the 
main purposes of this article is to report the first results for A hypernuclei calculated with a 
version of the QMC model extended to flavor SU(3), which takes account of the quark structure 
of the bound A as well as the bound nucleons. Furthermore, we treat both the A hyperon sitting 
outside of the closed-shell nuclear core, and the nucleons in the core, fully self-consistently, in 
the relativistic mean field approach. Extensive results for A, £ and S hypernuclei together will 
be reported elsewhere 0. 

Within the Born-Oppenhcimer approximation, one can derive equations of motion for a 
A hypernucleus in the QMC in the same way as has been done for normal nuclei |l], 0. Such an 
approach may provide us with important information about the hyperon-nucleon interaction, 
the deep nuclear interior and a possible manifestation of the quark degrees of freedom via 
the Pauli principle at the quark level. As an example, the very weak spin-orbit interaction 
for A hypernuclei, which had been phenomenologically suggested by Bouyssy and Hiifner [0- 
was first explained by Brockman and Weise [11| in a relativistic Hartree model, and directly 
confirmed later by experiment |12|. However, a very strong SU(3) breaking effect was required 
to achieve this small spin-orbit force. An explanation in terms of quark and gluon dynamics 
was made by Pirner |T3|. Alternatively, Noble [|T4| showed that this smallness of the spin-orbit 
force could be also realized without any large breaking of SU(3) symmetry, if an toAA tensor 
coupling was introduced, analogous to the anomalous magnetic moment of the A. However, 
Dover and Gal ]15| questioned whether the tensor coupling of the to meson to the A could be 
related to the anomalous magnetic moment, because the spin of A is entirely carried by the s 
quark in a naive SU(6) valence quark model, and this s quark couples exclusively to the </> meson 
according to the OZI rule. Later, Jennings 0] pointed out that within Dirac phenomenology 
the tensor coupling of the to meson to the A could be introduced in such a way as to guarantee 
that the direct to coupling to the spin of the A was zero - as one would expect from a simple 
quark model. The resulting spin-orbit force agrees with the result obtained in the present work 
on the basis of an explicit treatment of the quark structure of the A moving in vector and scalar 
fields that vary in space. Many other studies of the properties of hypernuclei have been made 
using relativistic, mean held theory (0-[0|. There has also been a great deal of experimental 


work 0, 0, 


As a second example, it has also been discovered that there is an overbinding problem in 
the light A hypernuclei, and the existence of a repulsive core or the necessity of a repulsive three- 
body force have been suggested to overcome the problem |I3| . The origin of this overbinding, 


which could not be explained easily in terms of traditional nuclear physics, was ascribed to 
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the Pauli principle at the quark level by Hungerford and Biedenharn [[D]. Investigations of 
this repulsive core in the A-nucleus system have been made by Takeuchi and Shimizu, and 
others |35], |36|, based on a nonrelativistic quark model. 

In the light of these earlier investigations, it seems appropriate to investigate the (heav¬ 
ier) hypernuclear systems quantitatively using a microscopic model based on quark degrees 
of freedom. For this purpose, the QMC model (which is built explicitly on quark degrees of 
freedom) seems ideally suited, because it has already been shown to describe the properties of 
finite nuclei quantitatively. 

Using the Born-Oppenhcimer approximation one can derive mean-held equations of mo¬ 
tion for a hypernucleus in which the quasi-particles moving in single-particle orbits are three- 
quark clusters with the quantum numbers of a hyperon or a nucleon. One can then construct 
a relativistic Lagrangian density at the hadronic level |[T]. |2|, similar to that obtained in QHD, 
which produces the same equations of motion when expanded to the same order in v/c : 


£qmc ~ ^n( r) 


t n p 

'3 , M i _iV'i 


H-d- M* N (a) - (+ g p —b{r) + -(1 + r 3 )A(r) )y 0 


07v(f) 


+ ^ Y if) [ij-d- My (a) - ( g u w(r) + g p I 3 6(f) + eQ Y A(r) boj 0y(f) 
~ ^[(Vcr(f)) 2 + m 2 a(r) 2 ] + ±[(Va;(f)) 2 + m^(f) 2 ] 

+ ^[(V6(f)) 2 + m 2 6(F) 2 ] + i( VA(r )) 2 , 


( 1 ) 


where V ; n (f) C0r(f)) an d 6(f) are respectively the nucleon (hyperon) and the p meson (the 
time component in the third direction of isospin) fields, while m CT , m u and m p are the masses 
of the er, t o and p mesons. g u and g p are the cn-N and p-N coupling constants which are related 
to the corresponding (u,d)-quark-u;, g^, and (u,d)-quark-p, g q p , coupling constants as g u = 3 g q 
and g p = g q p [|, |. 

In an approximation where the cr, u; and p mesons couple only to the u and d quarks 
(ideal mixing of the a; and 0 mesons, and the OZI rule are assumed), the coupling constants 
in the hyperon sector are obtained as g 4 = (n 0 /3)g OJ , and g p = g p = g q p) with n 0 being the 
number of the valence u and d quarks in the hyperon Y. and Qy are the third component 
of the hyperon isospin and its charge, respectively. The held dependent cr-N and cr-Y coupling 
strengths predicted by the QMC model, g a (a) and g\ (cr), are defined by, 


M* n {(J ) = M n - g a (a)cr(r), (2) 

My(a) = My - g^(a)a(r), (3) 


where Mj y ( My ) is the free nucleon (hyperon) mass. Note that the dependence of these coupling 
constants on the applied scalar held must be calculated self-consistently within the quark model. 
Hence, unlike QHD, even though g^/g a may be | in free space, this will not necessarily be the 
case in nuclear matter. More explicit expressions for g^ (a) and g a (cr) will be given later. From 
the Lagrangian density in Eq. (jl]), one gets a set of equations of motion for the hypernuclear 
system, 




[i'y ■ d - M* 


N 


b) - (gMf) + gp~^~b(r) + e -{l + t 3 )A(j) ) 7 o]0jv(f) = o, 


[b ’ & - M*(a) - (g%u(r) + g p I 3 b(r) + eQ Y A(r) )yo ]Vr(f) = 0, 

, __o 2 \ t dM^(c r) dMyia ) Y 

(-V 2 + ml)a(r) = -[ £ ’ }p s {f) - [ ^ ’ }p Y s (f), 


(4) 

(5) 


da 


2 


da 











= 9 aC N (cr)p s (f) + g a C Y (a)p s 0 r.), 
V? + ml)(j(r) = g u p B (f) + glp Y B {r), 

-V r 2 + m 2 p )b(r) = ^p 3 (r) + g*I% p\{r), 

-V 2 r )A(r) = ep p (r) + eQ Y p Y B (r), 


( 6 ) 

(7) 

( 8 ) 
(9) 


where, p s (r) (p Y (f)), pnij) (Pb{t)), / 23 (f) and p p (f) are the scalar, baryon, third component of 
isovector, and proton densities at the position r in the hypernucleus [jp], |8j. On the right hand 
side of Eq. a new, and characteristic feature of QMC beyond QHD [19|, |29[ appears, 


9M*(a) 


= 9 trC n (<t) and 


dM*{a) 


namely, da - — da 

Y1 a = 0). The effective mass for the hyperon Y is defined by, 


= 9%Cy((t), where g a = g a (a = 0) and g Y = 


g 


BMP 


Y 


cr 


d 


da 


= -n 0 g q a I dy ^ u Ay)^ u Ay) = ~ n o 9 %S Y {c) = - g Y (a)a 


I bag 


da L' 


( 10 ) 


with the MIT bag model quantities |1], [2j, 


My (a) 
S Y (a) = 

n*(a) = 
C Y (a) = 


n 0 n'(a) + (3 - noJW-zy + 


R 


Y 


Q*(q)/2 + 7n* U)d (a)Rl-(Q*(a) - 1) 

0*(cr)(n*(c7) - l) + ml d (a)R* Y l2' 

\jx 2 + (R Y m* u 4 (a )) 2 , fl* = + (R^m s ) 2 , m* ud (a) = m u4 - g q a a(r ), 

: Sy(<7)/Sy(0), gl = n 09 ;Sy(0) = ^ 9 „Sy(0)/S K (0) = ^g„T y /K . (11) 


Here, z Y , B, x, x s and m U) d )S are the parameters for the sum of the c.m. and gluon fluctuation 
effects, bag pressure, lowest eigenvalues for the (u, d) and s quarks, respectively, and the 
corresponding current quark masses with m u = rrid■ Zn and B (z Y ) are fixed by fitting the 
nucleon (hyperon) mass in free space. 

The bag radii in-medium, R% Y , are obtained by the equilibrium condition 
dM^ Y (a) / dR]y :Y \ji N y =r* n Y — 0- The bag parameters calculated and chosen for the present 
study are, (z^, z\) = (3.295, 3.131), (Rn, Ra) — (0.800, 0.806) fin (in free space), B f 4 = 170 
MeV, (m u , m d , m s ) = (5, 5, 250) MeV. The parameters associated with the u and d quarks 
are those found in our previous investigations ||. The value for the mass of the s quark in 
the MIT bag was chosen to be 279 MeV, in order to reproduce the mass of the A in that 
model [0. However, the final results turn out to be insensitive to the values of many of these 
parameters [I], |2]. The value for V Y / N turned out to be almost unity for all hyperons, so we 
can use T y /n = 1 in practice 0. At the hadron level, the entire information on the quark 
dynamics is condensed in Cjv^a) of Eq. (!6]). Furthermore, when this Cn,y{ct) = 1, which 
corresponds to a structureless nucleon or hyperon, the equations of motion given by Eqs. (||)- 
(||) can be identihed with those derived from QHD [0 0, 0, except for the terms arising from 
the tensor coupling and the non-linear scalar held interaction introduced beyond naive QHD. 
The parameters at the hadron level, which are already fixed by the study of infinite nuclear 
matter and finite nuclei are as follows: m^ = 783 MeV, m p = 770 MeV, rn a = 418 MeV, 
e 2 /47r = 1/137.036, g 2 / 47t = 3.12, g 2 /47r = 5.31 and g 2 / 47r = 6.93. 

For practical calculations, it has been found that the scalar densities C Y:Y (a) can be 
parametrized as a linear form in the a held |l|, [2|, || , 


C NjY (a) = 1 - a NjY x (g a a(r)), 


( 12 ) 
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where g a a = (3^S , Ar(0))a in MeV, and the value for the A hyperon is, a A = 9.25 x 10~ 4 
(MeV -1 ) for m u = 5 MeV, m s = 250 MeV, Rn = 0.8 fm ~ J? A . This parametrization works 
very well up to about three times normal nuclear density, pb — 3p 0 , with p 0 — 0.15 fm -3 . Using 
this parametrization, one can write down the explicit expression for the effective mass of the 
hyperon Y : 


My(a) = M Y - gl (cr)o-(r) ~ M y - jg a 

The origin of the spin orbit force for a composite nucleon moving through scalar and 
vector fields which vary with position was explained in detail in Ref. |[]J - c.f. sect. 3.2. The 
situation for the A is different in that, in an SU(6) quark model, the u and d quarks are coupled 
to spin zero, so that the spin of the A is carried by the s quark. As the a-meson is viewed here 
as a convenient parametrization of two-pion-exchange and the a; and p are non-strange, it seems 
reasonable to assume that the a, a; and p mesons couple only to the u and d quarks. The direct 
contributions to the spin-orbit interaction from these mesons (derived in sect. 3 of Ref.|l[) then 
vanish due to the flavor-spin structure. Thus, the spin-orbit interaction, Vg Q ( r)l ■ s, arises 
entirely from Thomas precession: 


1 - y (fU^)) 


<7 r . 


(13) 


U 


s.o.(r,)l-s = ~v A x 


dv A 
dt 


■ s = 


2 M; 


*2 i 


ryr 


y Wl{r) + gtu{r)] j l ■ s, (14) 


where, va = p A /M£, is the velocity of the A in the rest frame of the A hypernucleus, and the 
acceleration, dv\/dt, is obtained from the Hamilton equations of motion applied to the leading 
order Hamiltonian of the QMC model |l| . Because the contributions from the effective mass of 
the A, M%(r), and the vector potential, g^uj (r), are approximately equal and opposite we quite 
naturally expect a very small spin-orbit interaction for the A in the hypernucleus. Although 
the spin-orbit splittings for the nucleon calculated in QMC are already somewhat smaller than 
those calculated in QHD , we can expect even much smaller spin-orbit splittings for the A in 
QMC, by comparing the expressions obtained for the spin-orbit potentials between the nucleon 
and A [T]]. (See also Fig. [|.) In order to include the spin-orbit potential of Eq. (|I4|) correctly, 
we added perturbatively the correction due to the vector potential, — ,, . (±g*u(rj) f- s, 

to the single-particle energies obtained with the Dirac equation, Eq. (|5]). This is necessary 
because the Dirac equation corresponding to Eq. O leads to a spin-orbit force which does not 
correspond to the underlying quark model, namely: 

V s \(ry- s = (4[MJ(f) - S >(r)]) f- S. (15) 

This correction to the spin-orbit force, which appears naturally in the QMC model, may also 
be modelled at the hadronic level of the Dirac equation by adding a tensor interaction. 

We are now in a position to discuss the results. As already mentioned, the calculations 
are fully self-consistent for all fields appearing in Eqs. (fljMjQ 1 ). In principle, the existence of 
the A outside of the nuclear core breaks spherical symmetry, and one should include this in a 
truly rigorous treatment. We neglect this effect, since it is expected to be of little importance 
for spectroscopic calculations |37|, [38|| . However, we do include the response of the nuclear core 
arising from the self-consistent calculation, which is significant for a description of the baryon 
currents and magnetic moments, and a pure relativistic effect |38], |39]. We will discuss later 
the self-consistency effects between the nuclear core and the A. 
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In Fig. |T|, we show the effective masses of the nucleon and A as well as the baryon densities 
calculated for (a): \ 7 0, (b): ^Ca and (c): ^ 09 Pb. The results are for the lsi/ 2 A state, where 
effects of the A on the whole system are expected to be the largest. The effective masses in 
the hypernuclei, \ 7 0, Ca and ^ 09 Pb behave in a similar manner as the distance r from the 
center of each nucleus increases (the baryon density decreases). The decrease of the baryon 
densities around the center in \'0 in Fig. |l](a), which one might think a shortcoming of the 
present treatment, is also expected in nonrela.tivistic potential model calculation p0| . 


The calculated A single-particle energies for the closed-shell nuclear core plus one A con¬ 
figuration are listed in Table [l|. with the experimental data |2t], Q. At a glance, one can 


easily see that the spin-orbit splittings in the present calculation are very small for all hyper- 
nuclei. These small spin-orbit splittings, tend to be smaller as the baryon density increases, 
or the atomic number increases. This smallness of the spin-orbit splitting is a very promising 
achievement in the present approach. Concerning the single-particle energy levels, although 
a direct comparison with the data is not precise due to the different configurations, the cal¬ 
culated results seem to lead to overestimates. In order to make an estimate of the difference 
due to this different configuration, we calculated the following quantity. By removing one lp 3 / 2 
neutron in 16 0, and putting a A as experimentally observed, we calculated in the same way as 
for ]\ 7 0-which means the nuclear core was still treated as spherical. In this case the calculated 
energy for the lsi /2 A is -19.9 MeV, to be compared with the value -20.5 MeV of the present 
calculation. For larger atomic numbers, the difference is smaller. (See also Table [I].) We also 
attempted the calculation using the scaled coupling constant, 0.93 x g^(a = 0), which repro¬ 
duces the empirical single-particle energy for the lsi /2 in ACa, -20.0 MeV |29|. The results 
obtained using this scaled coupling constant are also listed in Table 1, denoted by ^ 4 Ca* and 
| 09 Pb*. Then, one can easily understand that the QMC model does not require a large SU(6) 
(SU(3)) breaking effect in this respect (7 %) to reproduce the empirical single-particle energies. 
Furthermore, the model still achieves the very small spin-orbit interaction for A in hypernu¬ 
clei, without introducing the tensor coupling of the u meson, which may be contrasted with 
the QHD type calculations |^2|. We note that the spin-orbit splittings for the A (or hyperon) 
single-particle energies in hypernuclei are not well determined by the experiments f28|-||32||. 

In Table we list the calculated binding energy per baryon, —E/A, rms charge radius, 
r c h , and rms radii of the A and the neutron and proton distributions (r\, r n and r p , respectively), 
for the lsi /2 and lp 3 / 2 A configurations. The rms charge radius is calculated by convolution 
with a proton form factor J2|. For comparison, we also give these quantities without a A-i.e., for 
normal finite nuclei. The differences in values for finite nuclei and hypernuclei listed in Table 
reflect the effects of the A, through the self-consistency procedure. One can easily see that the 
effects of the A become weaker as the atomic number becomes larger, and the A binding energy 
becomes smaller. 

Regarding the effects of the A on the core nucleons, we show also in Fig. |2| the core-nucleon 
single-particle energies for 40 Ca, and ^Ca for a lsi /2 A state. It is interesting to see that the 
effects of the A work differently on the protons and neutrons in ^ 4 Ca. The existence of the 
A makes the scalar and baryon densities larger, and the scalar and vector potentials become 
stronger. As a consequence, the relative strength between the scalar, vector and Coulomb 
potentials change in ^Ca, compared to those in 40 Ca. From the calculated results for the 
proton rms radius, r p , and nucleon energy shift shown in Fig. |2|, we can conclude that the 
Coulomb potential increases more than the gain in the difference between the scalar and vector 
potentials. We can also say that the scalar potential increases more than the vector potential, 
by observing the neutron single-particle energies and the rms radius of the neutron, r n . (See 
also Table §.) 
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Table 1: A single-particle energies 


in MeV) for closed-shell nuclear core plus one A config- 
are calculated using the scaled coupling constant 
0.93 x g£(cr = 0), which reproduces the empirical single particle-energy for the lsi /2 in ^Ca, 


uration. The values for ^Ca* 


and ^ 09 Pb* 


-20.0 MeV ^9f. For reference, we list the experimental data for \ 6 0 and ^°Ca of Ref. p9f de¬ 


noted by a, and for 89 Y and | 08 Pb of Ref. |32j denoted by b , respectively. Spin-orbit splittings 
are not well determined by the experiments. (* fit) 
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Table 2: Binding energy per baryon, — E/A (in MeV), rms charge radius, r c h , and rms radii of 
the A, r\, neutron, r n , and proton, r p (in fm). (* fit) 



A state 

-E/A 

f ch 

T A 

T n 

r p 

17 n 

A u 

lSl/2 

6.75 

2.85 

2.22 

2.58 

2.73 

17 o 

A u 

1P3/2 

6.09 

2.82 

2.95 

2.61 

2.70 

16 0 

5.84 

2.79 


2.64 

2.67 

a Ca 

lSl/2 

7.83 

3.52 

2.55 

3.30 

3.42 

£Ca 

1P3/2 

7.57 

3.51 

3.17 

3.32 

3.41 

40 Ca 

7.36 

3.48* 


3.33 

3.38 

! y Ca 

lSl/2 

7.75 

3.54 

2.61 

3.63 

3.45 

A 9 Ca 

lP3/2 

7.53 

3.54 

3.24 

3.64 

3.44 

48 Ca 

7.27 

3.52 


3.66 

3.42 

yl Zr 

lSl/2 

8.04 

4.29 

3.02 

4.29 

4.21 

91 Zr 

l/R/2 

7.98 

4.28 

3.66 

4.30 

4.21 

90 Zr 

7.79 

4.27 


4.31 

4.19 

7 J Pb 

lSl/2 

7.39 

5.50 

3.79 

5.67 

5.43 

20 9 pb 

1P3/2 

7.37 

5.49 

4.50 

5.67 

5.43 

208 pb 

7.25 

5.49 


5.68 

5.43 


Finally, we show the scalar and vector potential strength for \ 7 0, Ca and ^ 09 Pb in Fig. [|. 
The difference between the scalar and vector potentials near the center of the hypernuclei is 
typically ~ 35 — 40 MeV. This value is slightly (~ 5 — 10 MeV) larger than that calculated by 
Ma et al. |26|, which seems to be the specific origin of the overbinding in the present calculation. 
Although our calculations are based on quark degrees of freedom, in the effective hadronic model 
the whole dynamics of the quarks and gluons are absorbed into the parameters and coupling 
constants appearing at the hadron level. In this sense, the effect of the Pauli principle at the 
quark level |34|, [Pj, is not included between the quarks in the A and the nucleons. It seems 


necessary to incorporate the Pauli principle at the quark level, or some equivalent effect which 
would produce a repulsive core in hypernuclear systems, if one is to reproduce the experimental 
binding energies in the QMC model. 

In summary, we have reported the first results for A hypernuclei calculated with the QMC 
model extended to the flavor SU(3). The very small spin-orbit force for A in hypernuclei was 
achieved. This is a very promising feature and a direct consequence of the SU(6) quark model 
wave function and the quark structure of the A in the QMC model. However, the calculated 
single-particle energies for the A tend to be overestimated in comparison with the experimental 
data. In order to overcome this overbinding problem, it may be necessary to introduce a 
repulsive core due to the Pauli principle at the quark level, or some equivalent effect. This is a 
challenging problem for future work. 


We would like to thank R. Brockmann for helpful discussions concerning the calculation, 
and P.A.M. Guichon for useful comments. This work was supported by the Australian Research 
Council, and A.W.T and K.S. acknowledge support from the Japan Society for the promotion 
of Science. 
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Figure 1: (c) 


Figure 1: 

Calculated baryon densities, p#, and effective masses of the nucleon denoted by, N, and the A 
hyperon denoted by, A, in hypernuclei for (a): \ 7 0, (b): ^Ca and (c): ^ 09 Pb. All cases 

are for the lsi/ 2 A state. 
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Figure 2: Nucleon single-particle energies for 40 Ca, and ^ Ca for the lsi/ 2 A state 


12 




0 2 4 6 8 10 

r (fm) 


Figure 3: Calculated scalar and vector potential strength for ^'0, ^Ca and ^ 09 Pb. 
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